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Thomson scattering in high-intensity chirped laser pulses 
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We consider the Thomson scattering of an electron in an ultra-intense chirped laser pulse. It is found that the 
introduction of a negative chirp means the electron enters a high frequency region of the field while it still has a 
large proportion of its original energy. This results in a significant enhancement of the energy and intensity of 
the emitted radiation as compared to the case without chirping. 


I. INTRODUCTION 

In recent decades, since the discovery of chirped pulse am¬ 
plification [1], the powers and intensities of laser facilities 
around the globe have been exponentially increasing [2]. The 
current record of 2 x 10" W cm“ 2 was set in 2008 [3] and it is 
expected that this will be routinely surpassed in the near future 
at new facilities such as the Vulcan 20 PW upgrade [4], the Ex¬ 
treme Light Infrastructure (ELI) Facility [5] and the XCELS 
project [6]. The widespread availability of the current tech¬ 
nology has driven a large field of research in nonlinear Thom¬ 
son and Compton scattering, with a view to producing high- 
energy, tuneable y-ray beams. These sources have important 
applications, both for fundamental research [7] and for more 
practical applications, such as cancer radiotherapy [8] and the 
radiography of dense objects [9]. Recent experiments [10, 11] 
have been pushing the limits of peak energies and brilliances. 
However, working with ever-higher laser intensities we will 
soon enter a regime where radiation reaction (and, ultimately, 
QED effects [12, 13]) will start to come into play. 

Radiation reaction (RR) occurs when particles are acceler¬ 
ated so strongly by the laser field that their resulting radiation 
emissions cause significant energy losses. The result is a fric¬ 
tional effect which can significantly impact on the particle dy¬ 
namics, causing the particles to slow and reducing their energy 
as they reach the peak field [14, 15] (see also [16, 17]). As a 
consequence the resulting emission spectrum will be reduced 
in both energy and intensity (see, e.g. [18, 19]). 

Recently a number of articles have considered the effects 
of pulse chirping in laser-matter interactions. This has been 
in the context of ion acceleration [20], and in the Thomson 
scattering of relativistic electrons in moderately intense laser 
pulses [21, 22], In this paper we show how the introduction of 
a chirp into a very intense laser pulse can help to mitigate the 
reduction in energy of the emitted radiation by allowing the 
electrons to probe deeper into the laser focus before becoming 
susceptible to RR. The result is a higher electron energy in the 
peak field, enabling a significant increase in the energies and 
brillances of the Thomson radiation. 
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II. THEORY 

We consider the interaction of a relativistic electron with 
a counter-propagating chirped Gaussian laser pulse of base 
frequency ©o and FWHM duration To- We take the direction 
of the laser propagation to be along the z axis, such that the 
laser wave vector is k = ©oz/c. Then we normalize space and 
time with respect to the wave vector and the base frequency, 
respectively (x -A kx and t -A ©of). Defining a chirped pulse in 
the same manner as a number of recent works (e.g. Refs. [20- 
22]), the vector potential can be written as 

A = A 0 exp [5cos(?7+/(? 7 ))e. r 

+ \/l -<5 2 sin(77+/(t7))e v ], (1) 

where Ao is the field amplitude, £ = 41n(2), T) = t—z + <j)o, 0o 
is a phase constant, and /(rj) is the chirp function 1 . The factor 
8 determines the polarization, and is set to 1 (1/ y/2) for linear 
(circular) polarization. The vector potential is then normal¬ 
ized as A —► eA/me, and hence the electric field can be nor¬ 
malized as E —>• cE/mO)()C, where A and E are the vector po¬ 
tential and the corresponding electric field in SI units. We then 
define the dimensionless intensity parameter ao = eEo/mcOQC 
in the usual manor, in terms of the peak fields. 

A frequency chirped pulse is one where the frequency © 
changes with time. This can be achieved by the use of a set of 
grating pairs [23, 24], and it is now fairly trivial to introduce 
chirps such that the laser frequency changes by a few percent 
over the pulse duration. In our expressions the chirping of the 
laser pulse is defined by the function /(rj) = hi] 1 , where b is 
the chirping constant. The value of b must be kept small since 
the since the chirp is restricted by the bandwidth, and in turn 
the length, of the initial generating pulse. Setting b = 0 will 
correspond to the unchirped laser pulse. It is important to note 
that the introduction of a chirp into the laser pulse will change 
the pulse energy (~ |E| 2 ). In order to compensate for this we 
reduce the pulse duration accordingly. 


1 Note that introducing the chirp into the vector potential will result in dif¬ 
ferent electromagnetic field components, and therefore different particle 
dynamics, to those obtained by introducing the chirp directly into the E 
and B fields. 
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Having described the pulse, we now turn our attention to the 
motion of the colliding electrons. Ordinarily the particle mo¬ 
tion would be governed by the Lorentz equation, but in cases 
where the acceleration is strong the emission of radiation can 
lead to a significant reduction in the particle’s energy and mo¬ 
mentum. The effect of this ‘radiation reaction’ on the particle 
dynamics can be included by adding a correctional term to the 
Lorentz force equation. However, determining what this cor¬ 
rection should be is surprisingly non-trivial. Despite having 
been studied for over 100 years, it remains one of the most 
fundamental problems in electrodynamics. A common start¬ 
ing point is to solve the coupled Lorentz and Maxwell’s equa¬ 
tions for the system. This results in the infamous Lorentz- 
Abraham-Dirac equation [25-27], which suffers from notori¬ 
ous defects such as pre-acceleration and (unphysical) runaway 
solutions. A common resolution to these problems is to adopt 
the perturbative approximation of Landau and Lifshitz [28]. 
Then the equation of motion is given by 


dp 

dt 


= *L 


f L T fr? 


( 2 ) 


where ^ = Ef vxli is the Lorentz force and the radiative 
correction term 


fR 


4 r e 
3^1 


Jy+v-VjE + vx f+v - V ] B 


(f L )xB + (vE)E-r[(f L ) 2 



(3) 


where r e = e 2 /me 2 is the classical electron radius and A = 
271c/C0 q is the base wavelength of the field. Equation (3) is 
valid when the radiative reaction force is much less than the 
Lorentz force in the instantaneous rest frame of the particle. 
We stress that there are a number of alternative equations in 
the literature (for an overview see [29, 30]) and it is still an 
open problem as to which is the correct formulation. How¬ 
ever, it has recently been shown that the Landau-Lifshitz equa¬ 
tion, along with some of the others, is consistent with quantum 
electrodynamics to the order of the fine-structure constant a 
[31, 32], Also, we note that the first term (derivative term) 
of Eq. (3) is significantly smaller than the other two, since 
it is only linear in the field strength whereas the other terms 
are quadratic. It is found that in almost all cases the contri¬ 
bution from this term is negligible and so we do not include 
it in our simulations. (Indeed, it can be shown that, in cases 
where classical RR is important, the derivative term is even 
smaller than the electron spin force and so one could argue 
that it should be neglected out of consistency [33].) 

Once we have calculated the particle trajectory in the pulse, 
the resulting radiation emissions can be obtained via a well- 
known classical formula. The energy radiated per unit solid 
angle per unit frequency is given by [34], 


d 2 1 

dto d£l 


r nx[(n-p)x$] cisHD(t)]dt 
J (1 — /3 - n ) 2 


(4) 


where n is a unit vector pointing from the particle’s position to 
the detector ( D ) located far away from the interaction, and fi 



FIG. 1. Plots showing the time evolution of the electron y-factor in 
chirped and unchirped pulses (chirping constant b = ±0.008). The 
electron has an initial yo = 900 and the laser is a circularly polarised 
pulse of peak intensity ag = 10 and 10 cycles duration FWHM. 


and f> are, respectively, the particle’s relativistic velocity and 
acceleration. In our dimensionless units, s = co/cog is taken to 
be the harmonic of fundamental frequency. Here we have nor¬ 
malized the intensity by the factor e 2 /(47T 2 c). All the quanti¬ 
ties in the above equations are evaluated at the retarded time 
so one can directly do the integration in some finite limit. 


III. RESULTS 


Since this study is concerned with high field intensities, it 
is instructive to first provide an estimate for when RR effects 
become important. Using just the Lorentz force to determine 
the motion, the radiated power P is given by Larmor’s formula 
in terms of the acceleration, 

2 mr e 'dcc 2 2 , , . 

P=- -= -r e mca$<4Y(l+P), (5) 

3 c 3 

Normalizing this by come 2 we obtain the energy loss per cycle 
in terms of the electron rest energy me 2 [19, 35] 


R = 


P 2 too 2 


tOgmc‘ 


= -r e — 057 ( 1 +J3). 


( 6 ) 


When this parameter reaches unity we are in the “radiation- 
dominated regime” [36], where the radiation damping effects 
are of the same magnitude as the Lorentz force. 

We will begin by studying the effects of chirping on the 
Thomson spectra of a high-energy electron in a moderately 
intense laser, where RR effects are not very important. To be 
specific, we will consider an electron with initial yo = 900 (~ 
460 MeV) brought into collision with a circularly polarized 
laser pulse having peak amplitude ag = 10 (which corresponds 
to 2.11 x 10 20 W cm 2 for an optical laser). The duration of 
the unchirped (b = 0 ) laser pulse is taken to be 10 cycles and 
is reduced when the chirp is introduced in order to maintain 
a constant energy. For these parameters we find R=0.0027, 
meaning that we are in a regime where RR effects are likely 
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FIG. 2. Plots showing the time evolution of the electron y-factor for 
two different chirping constants. The electron has an initial yo = 900 
and the laser is a circularly polarised pulse of peak intensity ao = 200 
and 10 cycles duration FWHM. 


FIG. 3. Plots showing the emission spectra, evaluated in the 
backscattering direction (0 = 180°) for the two chirp factors (a) and 
(b). The emission spectra without RR taken into account is also pre¬ 
sented (c). Parameters as for Fig. 2 . 


to be minimal. We first solve the Landau Lifshitz equation for 
these parameters and consider how the y-factor changes as the 
electron passes through the laser pulse. In Fig. 1(b) we com¬ 
pare the evolution of the y-factor of an unchirped pulse with 
those of a positive and negative chirp of b = ±0.008. It can 
be seen that the electron in the negatively chirped field loses 
energy more quickly than that in the unchirped field, and the 
electron in the positively chirped field less quickly. (However, 
the overall energy loss is not large, amounting to about 5% 
of the starting value.) In the left hand panel. Fig. 1(a), we 
show the emission spectra for the three cases. This has been 
calculated by inserting the particle trajectories into (4). It can 
be seen that the electrons in the two chirped pulses emit ra¬ 
diation of a higher energy, but lower amplitude than the elec¬ 
tron in the unchirped pulse. This is not surprising since these 
electrons pass through regions of higher frequency fields than 
exist in the unchirped pulse. We will discuss this in more 
detail shortly. Note that the spectra for the two chirped cases 
are roughly similar. Observe also that the emission spectra are 
much cleaner for the electrons in the chirped pulses than in the 
unchirped pulse. This is due to the changing frequency caus¬ 
ing an overlap of the contributing harmonics and is discussed 
in Ref. [22], 

Having examined this preliminary example, we now move 
on to consider a case where RR effects do become important. 
For the rest of this study we will work with a circularly polar¬ 
ized laser pulse having peak amplitude oq = 200 (which cor¬ 
responds to 8.56 x 10 22 W cm~ 2 for an optical laser), which 
is just slightly beyond the current state of the art. The counter 
propagating electron will remain at yo = 900 (~ 460 MeV). 
For these parameters R = 1.06, placing us within the regime 
where RR effects can be expected to dominate. 

We consider once again how the y-factor changes as the 
electron passes through the laser pulse. In Fig. 2 we show 
the evolution of y for two different chirping constants, b = 
±0.002 and b = ±0.008. It can be observed, just as in the pre¬ 
vious example, that the energy of an electron colliding with a 


negatively chirped laser pulse (i.e. one where the high fre¬ 
quency components hit the particle first) falls more rapidly by 
virtue of RR as compared to both its positively (low frequency 
components first) or unchirped pulse counterparts. We note 
that the simple analysis used in eq. (6) is not sufficient to ex¬ 
plain this behaviour. There we were implicitly assuming that 
the chirped frequency ft) and y-factor can be approximated as 
constant over a laser cycle, which gives R ~ r e ye 2 |E| 2 / (Om 2 c 3 
for y 1. Thus, according to this simple model, increasing 
the pulse frequency will result in a smaller energy loss per 
cycle. However, the electron will pass through a correspond¬ 
ingly larger number of cycles during a given time period and 
thus the total energy loss over this time period will be approx¬ 
imately the same. Of course, this is only a crude estimate and 
we find that, when we calculate the particle motion numeri¬ 
cally, the higher the frequency components at the front of the 
pulse, the faster the particle loses energy. 

In Fig. 3 we consider the emission spectra, calculated by 



uj/ui 0 


FIG. 4. Angular distribution of the radiation spectra for b = 0 (a), 
b = —0.008 (b) and b = 0.008 (c). Parameters as for Fig. 2. 
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FIG. 5. Evolution of y along with laser pulse profile for b = —0.02 
(a) and b = 0.02 (b) is presented. Note that the chirp is unrealisti¬ 
cally large in this case and has been chosen to aid illustration. Other 
parameters the same as for Fig. 2. 


FIG. 6. Plots showing the time evolution of the quantum efficiency 
parameter, for two different chirping constants. Parameters are the 
same as for Fig. 2. 


inserting the particle trajectories into (4). At high intensities 
the total emission spectrum will be comprised of a sum of 
harmonics, corresponding to multiples of the laser frequency. 
In the unchirped case of constant frequency the properties of 
the spectrum are well understood. From conservation of mo¬ 
mentum arguments it can be shown that the frequency ft)' of 
the nth harmonic in the backscattering direction (9 = 180°) is 
given by [37] 




4 y 2 na>o 

1 ± og 


(7) 


For large «o the spectrum begins to decay after the critical 
harmonic with number « cr j t ~ 3 «q/2 [38], which for our pa¬ 
rameters is « cr i t ~ 10 7 . 

In our case the situation is more complicated than that of 
a monochromatic plane wave, but nevertheless we are able to 
make some very rough analogies to help us understand the 
physics. (For a detailed discussion of how the monochromatic 
spectra relate to those of a pulsed field we refer the reader 
to Refs. [39 —41 ].) We can reasonably assume that the spec¬ 
trum will be dominated by the emissions from the region just 
to the front of the pulse peak, where both the amplitude and 
the y-factor are relatively large. Assuming, for the purposes 
of this heuristic argument, that the field can be approximated 
over a given cycle by a monochromatic field with a the same 
amplitude, we can consider the above expressions in terms of 
the local dimensionless intensity a = e\E\/come. In particu¬ 
lar, we see that the number of harmonics comprising the spec¬ 
trum « cr j t will then be dependent on the local frequency. Thus, 
in the case of the positive chirp, although the frequencies of 
each harmonic comprising the spectrum will be blue shifted 
by the larger y-factor, the total number of harmonics will be 
decreased due to the lower frequency, resulting in a total spec¬ 
trum that covers a smaller energy range. Similarly, in the case 
of the negative chirp, the lower y-factor will red-shift the har¬ 
monic frequencies, but the total spectrum will be composed of 
a larger number of harmonics due to the higher frequency of 


the field. Thus the total spectrum will span a wider range of 
energies, as can be seen in the plots. Finally, the faster rate of 
energy loss for the negatively chirped case results in a stronger 
peak signal in the spectrum. 

The above argument can be supported by returning to our 
initial example where RR effects were insignificant. In this 
case no matter whether the chirp was positive or negative, the 
electron passed through the high frequency portion of the field 
with an energy comparable to its initial energy, without hav¬ 
ing lost much energy to RR. Thus it did not matter if it saw 
the high frequency part of the pulse first or later. This is why 
the radiation spectra was more or less identical for both cases 
(Fig. 1). Additionally, we can also consider the current ex¬ 
ample without the effects of RR (i.e. by solving the Lorentz 
equation instead of the Landau Lifshitz equation to determine 
the particle dynamics). The resulting spectra are plotted in 
Fig. 3(c) where it can be seen that the positive and negative 
chirp parameters produce identical spectra. 

Finally, a more complete presentation of the emission spec¬ 
tra is given in Fig. 4, for the chirping constant b = ±0.008. It 
can be seen from this figure that the spectra mainly consist of 
back scattered radiation, i.e. along 9 = 180°. 

In order to illustrate the the effect of the chirp more clearly, 
the interaction involving a slightly higher chirping constant 
(b = ±0.02) is presented in Fig. 5. In panels (a) and (b) 
we show how the y-factor of the electron evolves as it passes 
through the laser field. It can be seen that in the case of neg¬ 
ative chirp (panel (a)) the electron has a much higher y-factor 
in the high frequency part of the field than it does for the case 
of positive chirp (panel (b)). This results in a higher energy 
emission spectrum for the negative chirp compared to the pos¬ 
itive chirp. (Note that the results presented in Fig. 5 are only 
for the purpose of demonstrating the chirp. The energy loss in 
the early stages of the interaction with the positively chirped 
pulse is due to the ‘kink’ in the tail of the field. This is an 
unphysical artefact caused by our chirp being too large in this 
illustrative example.) 

Since the intensity we have been considering is fairly high. 
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it is worthwhile investigating whether it is indeed valid to treat 
the system classically, or whether quantum effects should be 
taken into consideration. We can measure the importance 
of quantum effects by considering the dimensionless ‘quan¬ 
tum efficiency parameter’ [42], which with our normalizations 
can be written % = hcooy\J (E + v x B ) 1 2 3 4 5 6 7 8 9 10 /mc 2 ~ yE /Eg, where 
Eg = 1.3 x 10 16 V/cm is the QED ‘critical’ field (‘Sauter- 
Schwinger’ field [43—45]). This parameter describes the work 
done by the laser field on the particle over of a Compton wave¬ 
length. When x ~ 1 quantum effects will start to dominate, 
with processes such as vacuum pair production occurring. Not 
only is the calculation of this parameter important in determin¬ 
ing the validity of our modelling, it is also of interest to see if 
the introduction of a chirp can significantly alter its value. We 
have already seen how the negative chirp results in an increase 
in radiated energy, and so it is worth investigating whether the 
same mechanism will increase the importance of QED effects. 
In Fig. 6 we plot the evolution of the % parameter for the cases 
we have been considering. It can be seen that % reaches a peak 
value of roughly 0.2 meaning that we are on the threshold of 
where QED effects are likely to be detectable, but not signif¬ 
icant [46], Thus we are justified in performing our analysis 
classically. We also see that, while the negative chirp does re¬ 
sult in a slight increase in j, the effect of chirping on QED 
effects is fairly minimal. This means that chirping is unlikely 
to prove a useful tool for probing intensity effects in strong 
field QED. 


IV. CONCLUDING REMARKS 

In this article we have investigated the dynamics, and result¬ 
ing Thomson spectra, of an electron in an intense chirped laser 
pulse. Because of RR effects the electron loses energy quickly 
upon entering the pulse. This means that by the time the elec¬ 
tron reaches the most intense part of the field it has a much 
lower energy than it began with. The result is a reduction 
in energy and brilliance of the emitted Thomson radiation. By 
introducing a small, negative chirp into the laser pulse we have 
shown that it is possible to have the electron enter the region 
of the pulse where the emissions will be strongest while it still 
has a large proportion of its initial energy. For the modest 
chirp parameters that we have considered this can result in a 
more than doubling of both the maximum frequency and am¬ 
plitude of the radiation spectrum a compared to the case of 
an unchirped field. This is of great importance in the context 
of Thomson/Compton scattering experiments using the next 
generation of ultra-intense laser sources. 
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